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1 , Abstract 

. The generation of an undular bore in the vicinity of a wave-breaking point is con- 

I sidered for the integrable Kaup-Boussinesq shallow water system. In the framework 

' of the Whitham modulation theory, an analytic solution of the Gurevich-Pitaevskii 

• type of problem for a generic "cubic" breaking regime is obtained using a generalized 

CN ! hodograph transform, and a further reduction to a linear Euler-Poisson equation. The 

I motion of the undular bore edges is investigated in detail. 

o : 

■ 1 Introduction 

^ i It is well known that the solutions to an initial value problem for the non-dispersive shallow 
water equations may lead to wave-breaking after a finite time, when the first spatial deriva- 



X 



^ • fives blow up. After the wave breaking point, a formal solution becomes multi-valued and 



loses its physical meaning. The divergence of the spatial derivatives at the wave-breaking 
point suggests that dispersion effects described by terms with higher order spatial deriva- 
tives must be taken into account. Then these small dispersion effects lead to the onset of 
oscillations in the vicinity of the wave breaking point followed by the development of an 
undular bore, or in different terminology, a dissipationless shock wave. 

This physical picture has been put into mathematical form for waves described by the 
Korteweg-de Vries (KdV) equation by Gurevich and Pitaevskii (GP) in [T] (see also Whitham 
in In the GP formulation, the region of oscillations is presented as a slowly modulated 
periodic wave solution of the KdV equation. The parameters of the wave change little on 
typical wavelength/period scales which permits one to apply the Whitham modulation theory 
(j21l3])- The resulting Whitham equations describe the slow evolution of the parameters in 
the undular bore. In the original paper two typical problems were considered. One 
problem is concerned with the description of decay of an initial discontinuity for the KdV 



1 



equation and the exact x/t-similarity solution of this problem was constructed. Another 
problem corresponds to the (universal) initial stage of development of a bore when the 
solution of the dispersionless equation can be approximated locally by a properly chosen 
cubic curve. This problem was studied in numerically. An exact analytic solution to this 
problem was later obtained by Potemin |1] using Krichever's algebro-geometrical procedure 
for integration of the Whitham equations 0, [7j. Later, Potemin's solution was put into the 
general context of Tsarev's generalized hodograph transform [Hj in fOlElEl, where the 
hodograph equations were reduced to the classical Euler-Poisson equation. 

However, the KdV equation describes unidirectional propagation of nonlinear dispersive 
waves. An integrable bi-directional analog of the KdV equation was derived by Kaup in 
jH] using the Boussinesq approximation for shallow water waves Like the KdV equa- 
tion, the Kaup-Boussinesq (KB) system is completely integrable and therefore a powerful 
inverse scattering transform method can be applied to its investigation. In particular, the 
multi-phase periodic solutions of the KB system were found in [HI, the Whitham theory of 
modulations was applied in 13j to the problem of the decay of an initial discontinuity, and 
a quasiclassical description of soliton trains arising from a large initial pulse was developed 
in [H]. 

In this paper, we further extend the Gurevich-Pitaevskii theory to the case of bi-directional 
shallow water equations using the KB system, and construct an analytic solution to the 
Whitham-KB equations for the regime of generation of an undular bore in the vicinity of a 
breaking point. The obtained solution, along with its own significance in the representation 
of undular bores, will serve as an intermediate asymptotic in a more general formulation we 
are presently undertaking, in which small dissipation is taken into account. 



2 Periodic waves in the Kaup-Boussinesq system 

In dimensionless units (see, e.g. ^2]) the KB system can be written in the form 

ht + {hu)x + \uxxx = 0, , s 

Ut + UUx + K = 0, 

where h{x, t) denotes the height of the water surface above a horizontal bottom and u{x, t) is 
related to the horizontal velocity field (at the leading order it is the depth-averaged horizontal 
field). 

The KB system (0) is completely integrable and can be represented as the compatibility 
condition of two linear equations jHj 

ij^, = Aij, = -^B.,ij + Bij^ (2) 



with 



1 \' . f. I 



A=[\--u\ -h, B=-[\ + -u]. (3) 



Thus, the inverse scattering transform method can be applied for its investigation. In par- 
ticular, the periodic solution of ((H) can be obtained by the well-known finite-gap integration 
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method (see, e.g. ^3]) in the following way. Let V'+ and be two basis solutions of the 
second order spatial linear differential equation Then their product 

g = (4) 

satisfies the third order equation 

g^xx - "iAxg - 4:Agx = 0. (5) 
Upon multiplication by g, this equation can be integrated once to give 

^ggxx-\9l-Ag' = -P{X) (6) 

where the integration constant -P(A) can only depend on A. The time dependence of g{x,t) 
is determined by the equation 

gt = Bgx - Bxg. (7) 
This equation can readily be put in the form 



9) t \9 



(8) 



which can in turn be considered as a generating function of an infinite sequence of conser- 
vation laws. 

The periodic solutions of the system (|T)) are distinguished by the condition that -P(A) 
in (jni) be a polynomial in A. The one-phase periodic solution, which we are interested in, 
corresponds to the fourth degree polynomial 

4 

P(A) = JJ(A - A,) = A^ - siA^' + S2X^ - ssA + s^. (9) 

i=l 

Then we find from Eq. (jH)) that g{x,t) is the first-degree polynomial, 

g{x,t) = X-fi{x,t), (10) 

where ^{x,t) is connected with u{x,t) and h{x,t) by the relations 

u{x, t) = Si — 2fi{x, t), h{x, t) = — S2 — + si/i, (11) 

which in turn follow from a comparison of the coefficients of A* on both sides of Eq. (jH)). The 
spectral parameter A is arbitrary and on substitution of A = into Eq. (jH)) we obtain an 
equation for /x, 

while a similar substitution into Eq. ((Tj) gives 

f^t = -(/i + ^U)fix = -^Siflx- 
Hence, ii{x,t) as well as u{x,t) and h{x,t) depend only on the phase 

e = x~^sit, (12) 
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so that V =lsi = ^Y.^i (13) 

i=l 

is the phase velocity of the nonhnear wave, and ^{6) is determined by the equation 

fie = (14) 

For the fourth degree polynomial Q the solution of this equation is readily expressed in 
terms of elliptic functions. Let the zeros Aj, i = 1,2, 3, 4, of the polynomial -P(A) be real and 
ordered according to the rule 

Ai < As < A3 < A4. (15) 

Then the real variable n oscillates in the interval where the expression under the square root 
in (fT^ is positive, 

X2<ii< A3. (16) 
Consequently the solution of Eq. ()14|) with the initial condition /i(0) = A3 is given by 



A3(A4 - A2) - A4(A3 - A2)sn2 f V(A4 - A2)(A3 - Ai) m 
A4 - A2 - (A3 - A2)sn2 (v/(A4-A2)(A3-Ai)e, 1 



where 

(A3-A2)(A4-Ai) 



m 



(A4-A2)(A3-Ai) ^ ' 

is the modulus of the elliptic functions. An equivalent solution corresponding to the initial 
condition /i(0) = A2 is given by 



A2(A3 - Ai) - Ai(A3 - A2)sn2 V(A4 - A2)(A3 - Ai) ^, m 

m = . ; = (19) 

A3 - Ai - (A3 - A2)sn2 (V(A4-A2)(A3-Ai)0, mj 

Substitution of (fTTj) or (fT^ into (fTTj) gives expressions for u{9) and h{9) in the periodic 
nonlinear wave. Its wavelength is given by 

^3 du 2K(m) , , 

^ ^ ' (20) 



A. /PR V(A4-A2)(A3-Ai) 

K(m) being the complete elliptic integral of the first kind. 

The soliton limit (m = 1) is obtained either for Ai = A2 or for A3 = A4. For Ai = A2 
Eq. gives 

m = K {K-Xr){K-\^ ^21) 

A4 - Ai + (A3 - Ai)/cosh2[v/(A4 - Ai)(A3 - Ai) 9] 

and for A3 = A4 Eq. (fTIH) gives 

m = Ai + (22) 

A2 - Ai + (A4 - A2)/cosh2[^(A4 - A2)(A4 - Ai) 9] 
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Their substitution into yields the sohton solution of the KB system. 

In the opposite limit A2 = A3 (m = 0) both expressions (|17p and (jl9p reduce to 

= A2 = A3 . (23) 

Thus the limit A2 —>■ A3 yields sinusoidal waves. 



3 Whitham modulation equations for the Kaup-Boussinesq 
system 

The Whitham modulation equations describe the slow evolution of the parameters A,, i = 
1, 2, 3, 4, of a modulated nonlinear wave. They are 

^ + ^.(A)|^ = 0, ^ = 1,2, 3, 4, (24) 

where the Whitham velocities fi(A) can be expressed in the form 

v.{X)= (^l--^d?jV, 9, = -^, 1 = 1,2,3 A, (25) 

where the phase velocity V and the wavelength L are given correspondingly by p3|) and 
(pUjl . A simple calculation yields the explicit expressions [THj, 





(A4- 


-Ai)(A2-Ai)K(m) 


(A2 


(A3- 


)K(m) + (A4 - A2)E(m) 
-Ai)(A2-Ai)K(m) 


(A2 


-Ai: 
(A4- 


)K(m) - (A3 - Ai)E(m) 
-A3)(A3-A2)K(m) 


(A4 


-A3; 
(A4- 


)K(m) - (A4 - A2)E(m) 
-A3)(A4-Ai)K(m) 


(A4 


-A3; 


)K(m) + (A3 - Ai)E(m) 



V2 
V3 

where K(m) and E(m) are complete elliptic integrals of the first and second kind, respectively. 
In the limit Ai = \2 {m = 1) the Whitham velocities reduce to 

Vi=V2 = ]^'^\i, W3 = ^(3A3 + A4), = i(A3 + 3A4); (27) 
in the limit A3 = A4 (m = 1) they reduce to 

i^i = ^(3Ai + A2), i;2 = ^(Ai + 3A2), v^ = Vi = ]^Y.\i- (28) 
and in the limit A2 = A3 (m = 0) they reduce to 

-i = ^(3A, + A4), .2 = ^3 = ^E^- + '^t'/A!-'2A2^'^ ' -4 = ^(A. + 3A4). (29) 

Next we shall apply the Whitham theory to the description of the undular bore forming 
in the vicinity of a wave-breaking singularity. 
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4 The Gurevich-Pitaevskii problem for the KB system 



4.1 Wave breaking in the dispersionless limit 

In the dispersionless limit, the KB system reduces to well-known shallow water equations 

ht + (hu)^ = 0, Ut + uu^ + = 0, (30) 
which can be transformed to the diagonal form 

where X± = ± Vh (32) 

are the Riemann invariants of Eqs. ()3()|1 . 

Initial data are given by two functions A4.(x, 0) and A_(x, 0) determined by the initial 
distributions hQ{x) and Uo{x). The system i\'Sl\i has two families of characteristics in the {x, t) 
plane along which one of two Riemann invariants (either A+ or A_) is constant. The wave- 
breaking point corresponds to the moment when characteristics of one of the families begin 
to intersect, so that the corresponding Riemann invariant becomes a three- valued function 
in the physical plane. Let such an intersection occur for the characteristics transferring 
the values of A+. Then at the wave-breaking point the profile of A+ as a function of x 
has a vertical tangent line and, hence, in vicinity of this point it varies very fast, whereas 
the second Riemann invariant varies with x more slowly and may be considered here as a 
constant parameter: 

A_ = Ao = const. (33) 

Thus, in the vicinity of the breaking point we are dealing with a simple wave. The second 
equation in (jHH) is identically satisfied by Eq. (jHSl)- The first equation in (pTTjl then has the 
well-known solution 

^_1(3A^ + A_)t = /(A+), (34) 

where /(A+) is an inverse function to an initial profile A+(x,0). At the wave-breaking time, 
normalized here to be t = 0, the function x = f{X+) must have an inflexion point with a 
vertical tangent line. Then in the vicinity of this point /(A+) can be approximated by a 
cubic function, 

x-^(3A+ + A_)t = -C(A+-Ai)3. (35) 
The KB system is invariant with respect to Galilean transformation 

x' = X — Uot, t' = t, h = h\ u = u' + Uo, A± = A^ + (36) 

and scaling transformation 

x = ax', t = o?t\ h = h'/a?^ u = u/a, A± = A'^/a. (37) 
With the aid of these transformations Eq. ()35|) can be cast into the form 

x-^(3A+ + A_)t = -A^, A_ = Ao. (38) 
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Figure 1: Wave breaking of the water elevation in the dispersionless hmit; A_ is taken equal 
to -10. 

where we have omitted the prime superscripts for notational convenience. It corresponds to 
the wave breaking picture shown in Fig. 1. 

The actual solution of the KB system now consists of two parts. Following Gurevich and 
Pitaevskii [T], we suppose that the region of oscillations can be approximated by a modulated 
periodic solution of the KB system. Its evolution is determined by the Whitham equations 
()24|) and we have to find that solution which matches the solution (jHHjl at the end points of 
the oscillatory region. One may say that this oscillatory region (the undular bore) "replaces" 
a non-physical multi- valued region of the solution (jHHj) . One should emphasize, however, that 
the boundaries of the undular bore do not coincide with the boundaries of formal multi- valued 
solution. Outside these boundaries, the solution approaches the dispersionless solution (jHEI)- 

4.2 Undular bore solution 

We look for the solution of the Whitham equations ()24|) in the form 

X -v,{X)t = Wi{\), 2=1,2,3,4, (39) 

where Vi are the Whitham velocities (j26|) . Since we consider the breaking of the Riemann 
invariant A+ and A_ < A+, we take 

Ai = A_ = Ao = const. (40) 

Then the limiting formulas ()27|) - (P^ show that if we find Wi{\) such that 

wa = -A4 at A2 = A3 (m = 0), 
W2 = — A2 at A4 = A3 (m = 1), 

then Eqs. fl39|) determine the Riemann invariants in such a way that A4 = A-|. at the trailing 
edge x~{t) where m = 0, A2 = A_|_ at the leading edge x~^{t) where m = 1, and Ai = A_ = Aq 
everywhere. Thus, the plots of Riemann invariants A2, A3, A4 as functions of x are joined 
into continuous curve whose upper and lower branches match with the solution (j38|) of the 
dispersionless equations (see Fig. 2). In the region x~{t) < x < x~^{t) there are four Riemann 
invariants which determine the modulated periodic solution representing the undular bore. 
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Figure 2: Dependence of Riemann invariants A2, A3, A4 on x at fixed value of time t = 1 and 
with Ai = —10. The dashed hne shows the corresponding dependence of A+ for the solution 
of the KB equations in the dispersionless limit. 

At its trailing edge x x~{t) the amplitude of oscillations vanishes and at the leading edge 
X — i> x^{t) the periodic solution transforms into a soliton train. 

According to the generalized hodograph method [El EE], Eqs. ()39|1 satisfy the Whitham 
equations provided Wi{\) are the velocities of the flows 

^ + «;,(A)^=0, z = 1,2,3,4, (42) 

commuting with (j^ . i.e dtrXi = drtK- If "we represent Wi{\) in the form analogous to 
Eqs. (EH), 

L 



WiiX) = (^1 - —d^j W, z = 1, 2, 3, 4, (43) 

then the condition of commutativity of the flows ()24|1 and ()43|) reduces to the system of 
Euler-Poisson equations, exactly as happens in the KdV [10^ illj il2j and NLS [T7j cases, 

didjW- ^ (d,W-djW)=0, 1^3. (44) 



It is easy to check that this equation has a particular solution W = const/ a/ -P(A), -P(A) = 
n(A — Aj), which is sufficient for our purpose. We choose the normalization factor so that 
the coefficient before A~^ in the series expansion of W in powers of A~^ be equal to the phase 
velocity of the periodic wave Si/2 = V. Thus, we obtain the sequence of W^''^ defined by 
the generating function 

^ = = E ^ = 1 + ■ 1 + iM - 1^2) • ^ + iM - is,s, + iss) ■ ^ + • • • ■ 

(45) 

Next a sequence of velocities of the commuting flows is given by 

wi'\X) = (^1 - ^d.^ W^'\ z = l,2, 3, 4, (46) 
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where w) 



(1) 



Vi coincide with the Riemann velocities 



It is not difficult to find the 



limiting formulas analogous to (j?rj) - ()29p . In particular, we get at A3 = A4 (m 



(1) 



w 



(2) 



(1) 



w 



(3) 



A3=A4 



A3=A4 



A3=A4 



(2) 



(3) 



A3=A4 



A3=A4 



A3=A4 



2(^1 + 3A2), 

|(A? + 2AiA2 + 5A: 



2>i 



(47) 



— (5A? + 9A?A2 + 15AiA^ + 35A 
lb 



Now we take such the linear combination 

, (1) , (2) , (3) 

W2 = aQ + aiW2 + a2W2 + Ci^W2 

so that W2 satisfies the condition (PT|) . The coefficients 01,02,03,04 depend on the constant 
Riemann invariant Ai = A_ = Aq and their values found in this way yield the required 
solution of the Whitham equations: 



X — Vi 



(A)t = -f,wf\\) + |Ao«;r + iAo\(A) + ^Ao 
Ai = Ao = const. 



2,3,4; 



(48) 



These formulas define A2,A3,A4 implicitly as functions of x and t and give the solution of 
the Gurevich-Pitaevskii problem for the KB- Whitham system. It is interesting to note that, 
unlike the KdV case, this solution is not scale-invariant (we recall that in the counterpart 
KdV solution A^ = t-^'Hi{x/t^/^),i = 1, 2, 3). This happens due to the presence of the fourth 
Riemann invariant Ai in the Whitham equations which is constant (Aq) for the obtained 
solution and cannot be eliminated by simple Galilean transform. As a result, the solution 
(|48p does not possess the scaling invariance required for a generalized similarity behaviour 
of \j , j = 2, 3, 4. The only family of admissible such similarity {x/t) solutions is realized in 
the simplest case of the decay of an initial discontinuity studied in 



4.3 Laws of motion at the trailing and leading edges of the oscil- 
latory region 

Let us find the laws of motion at the leading and trailing edges of the undular bore. First 
we consider the leading edge A3 = A4 (m = 1), and define the small deviations A3 and A4 
from the value A3 = A4 = A4 : 

A3 = A+ + A'3, A4 = A+ + a; . (49) 
Then we seek the asymptotic expressions for formulas PH|) with i = 3,4 for small IA3I, IA4I: 
x~^ + x' — {v^ + v'^)t = + w'^, x~^ + x' — {v^ + v'^)t = + (50) 
where x' denotes the space coordinate reckoned from its limiting value x~^ and 

vt = vt = l{Xi + X2 + 2X4), (51) 
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- <^ a; In 



-(A2-A0A^ 
(A4-Ai)(A4-A2 



-a; In 



(A2 - Ai)A^ 



(A4-Ai)(A4-A2) 



Wo 



< = -TT7(5A^ + 6A^A4 + 
35 



iAaA^ + lGA^), 



(52) 
(53) 



Wo 



35 



(3A^ + 8A2A4 + 24A^)x 



-(A2 — Ai)A3 



a;, In 



(A4-Ai)(A4-A2) 



Alln 



(A2-Ai)A^ 
(A4-Ai)(A4-A2) 



(54) 



while Ai, A2, A4 denote here their hmiting values Xi, A^, A4 , correspondingly. Then subtrac- 
tion of one equation (jHn|) from the other yields at once the relationship 



^= ^(3A^ + 8A2A4 + 24A2). 



(55) 



On the other hand, the limiting formulas 



with 



5(Ai + 3A2) 



wt = -As 



and v^, given by and (f^Hj) give after subtraction the relationship 

t = 4(15A^ + 12A2A4 + 8A^). 
35 



(56) 



(57) 



Then equating of the right hand sides of (f^ and (fS7j) shows that at the leading edge we 
have 

3 



A 



--A+. 
4 ^ 



Substitution of this relation into or (j37j) gives the dependence of A2 on t: 



1/2 



(58) 



(59) 



At last, the formula — v^t = — (A^)^ yields the law of motion for the leading edge: 

xHt) = lXot+l^-t'^'. (60) 
In a similar way we can consider the trailing edge A3 = A2 = A^ [m = 0) where we define 

A3 = A2~ + a;, (61) 



A2 — A2 + A2, 



and Eqs. (j48j) with i = 2,3 reduce to 

X~ + x' — + V2)t = W2 + w'2, X~ + x' — {v^ + v'^)t = + W3, (62) 
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where 

_ _ A? + 4A1A2 - 8A^ - 2A1A4 + 4A2A4 + 
-""^ - 2(Ai - 2A2 + A4) ■ 



(63) 



A' + 3A' 

^2 = -^3 = -T77r^^r^TT(3^? - 8^i^2 + 8X1 + 2A1A4 - 8A2A4 + 3X1), (64) 

Z(^Ai — ZA2 + aa) 



= W3 =[128Xt - 64A^A4 - 16A^A^ - SAaA^ - 5A^ 

- 7Ai(16A^ - 8A^A4 - 2A2A2 - A^)]/(35(Ai - 2A2 + A4)) 



w' = -w' = ^ \ ^ , ,J - 384A^ + 384A^A4 - 80A^A^ - 16A2A^ - 9A^ 

' ' 70(Ai - 2A2 + A4)2 ^ ' ' ' ' ^ ' (66) 

+ 7A2(-24A2 + 8A2A4 + A^)]. 



(65) 



Then subtraction of one equation (jU^ from the other gives 

t = A[ _ 384A^ + 384A^A4 - 80A^A^ - 16A2A^ - 9A^ 

+ 7A^(-24A^ + 8A2A4 + A^)]/(3A2 - 8A1A2 + 8A^ + 2A1A4 - 8A2A4 + 3A^ 

Now the hmiting formulas 



(67) 



X — -{Xi + 3Xi)t = —Xl, X — v^t = —w^ (68) 

give 

^ _ 2 (8A2 - 7Ai)(8Ai + 4A2A4 + 3A1) - ISAj 

35 4A2 - 3Ai - A4 ■ ' 

Equating the right hand sides of ()67|1 and (jEHI), we find the relationship between the values 
of A2 and Aj at the trailing edge: 

21A?(A4 + 4A2) - 10Ai(20A^ + 2A2A4 - A^) + 16(8A^ - A^A4 - X2XI) + 9A^ = 0. (70) 

Given t and Ai = Aq = const, we find A2 = A2^ and A4 = A4 from (jHTIjl and (ffOjl and then the 
law of motion of the trailing edge follows from 

x- = ^(Ai + 3A4)t- Ai (71) 

It is worth noticing that Eq. (|HH|l coincides with the corresponding relation between the 
Riemann invariants at the leading edge in the solution of the Gurevich-Pitaevskii problem in 
the KdV equation case (see, e.g. J3]). A similar relation A^ = — Aj /4 between the Riemann 
invariants of the KdV theory follows from (fTOj) in the limit |Ai| = |Ao| ^ c)o. In the next 
approximation we obtain 



and similar expansion of Eq. ()69|) in powers of 1/A 



I7 



1,2 A , 10A4 



(-jAi 1 + -^), (73) 
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Figure 3: Riemann invariant A2 as a function of A4 defined implicitly by Eq. ()70p . The plot 
corresponds to a fixed value of Ai = —10. Dashed line shows the dependence according to 
asymptotic formula (fTij). 

yields with the same accuracy 

A4^V3t-^^, V3t<|Ai|. (74) 
7 Ai 

Dependence of A2 on A4 given by Eq. (fTOj) with fixed value of A 1 is illustrated in Fig. El 

Substitution of Eq. (|73|) into Eq. (|7ip yields an approximate expression for the law of 
motion of the trailing edge; 

x- = iAot-^t^/^ + ^|^, v^«|Ao|. (75) 

Thus, the analytic formulas ()48p for the solution of the Whitham equations allowed us to 
find the main characteristics of the dissipationless shock. With the use of Eqs. (j48|l we 
can find A2, A3, A4 as functions of x at given t, and their substitution into (fTTj) and (fTTjl 
yields the profiles of u{x) and h{x) in the undular bore. An example of such a profile of 
the water elevation is shown in Fig. |3J As we see, the non-physical solution obtained in 



h 



/ 


35 

N 30 

\ 

/ 

\ 25 




-10 -8 -6 ^ 


-\c4^__^ -2 









X 



Figure 4: Undular bore for the KB equations. The plot corresponds to the time t = \ and 
Ai = —10. Dashed line shows the solution in the dispersionless limit. 
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the dispersionless limit is replaced by an undular bore. Its end points move according to 
the laws found above, so that the oscillatory region expands with time and its width grows 
mainly as t^l'^. Small amplitude oscillations are generated at the trailing edge (actually 
these oscillations represent gravity waves propagating into an undisturbed smooth region), 
and they transform gradually into solitons at the leading edge. 

5 Conclusion 

In this article, the wave breaking (Gurevich-Pitaevskii) problem is solved for the shallow 
water waves described by the Kaup-Boussinesq system. This theory generalizes the KdV 
model to the case of a bi-directional wave propagation model, and so allows for the effects 
of wave interaction. As a result, the dispersionless theory includes two Riemann invariants 
and the breaking of the wave means breaking of one of these two invariants. Thus, the 
wave breaking picture depends on the value of an additional parameter: the value of the 
non-breaking Riemann invariant. Correspondingly, the Whitham theory for the KB system 
contains four Riemann invariants (instead of three in the KdV case) and, as a result, the 
solution of the Whitham equations for the Gurevich-Pitaevskii problem of the wave front 
breaking is parameterized by a constant value Aq. 

The analytic solution of the Whitham equations obtained here allows one to investigate 
the form of an undular bore as a function of the parameter Aq. In particular, the velocity of 
the trailing edge increases with decrease of |Ao| and in the limit |Ao| — *■ oo the present theory 
reduces to known results of the KdV model. 

Our obtained solution can also be viewed as an intermediate asymptotic in a more general 
problem of the description of frictional shallow water undular bores where small dissipation 
is taken into account (see [IHl [131201 for the results relevant to the KdV equation with weak 
dissipation). This problem will be the subject of a separate study. 
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